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The magnetic field (H) - temperature (T) phase diagram of a superconductor is significantly 
altered when domains are present in an underlying ferromagnet with perpendicular magnetic 
anisotropy. When the domains have a band-like shape, the critical temperature T c of the supercon- 
ductor in zero field is strongly reduced, and the slope of the upper critical field as a function of T 
is increased by a factor of 2.4 due to the inhomogeneous stray fields of the domains. Field compen- 
sation effects can cause an asymmetric phase boundary with respect to H when the ferromagnet 
contains bubble domains. For a very inhomogeneous domain structure, T c oc H 2 for low H and 
T c oc H for higher fields, indicating a dimensional crossover from a one- dimensional network-like to 
a two-dimensional behavior in the nucleation of superconductivity. 

PACS numbers: 74.25.Dw 74.25. Ha 74.76.Db 75.70.Kw 



I. INTRODUCTION 



In hybrid superconductor / ferromagnet (SC/FM) 
bilayers the FM modifies quite substantially the super- 
conducting properties of the SC layer. In particular, 
strong vortex pinning was reported recently for super- 
conducting films covering arrays of ferromagnetic dots 
with in-plan o 1 : 2 : 3 and out-of-plane magnetization, 4,5 
and for continuous SC/FM bilayers .SiguLiil Theoretical 
investigations showed that supercurrents and vortices 
can be induced in the SC by the stray field of the 
FM) 11 : 12 : 13 ! 14 : 15 : 16 and that the domain structure of soft 
FM's can be influenced by t he presence of the SC>i£ 
Furthermore. Ifladovic et alJ predicted the appearance of 
the so-called ir- phase state in SC/FM multilayers, where 
the phase of the superconducting order parameter tjj 
shifts by ir when crossing a ferromagnetic layer^S Re- 
cently the existence of the 7r-phase state was confirmed 
by observing sharp cusps in the temperature dependence 
of the critical current in SC/FM/SC junctions. 19 Earlier 
experiments were performed in order to find the 7r-phase 
state by measuring the predicted oscillatory dependence 
of the critical temperature T c of SC/FM multilayers 
on the FM layer thickness rf/ m i 20 i 21 ' 22 i 23 However, the 
results of these experiments were not conclusive, because 
the nonmonotonic T c (df m ) behavior could also appear 
due to the presence of magnetically "dead" layers at the 
SC/FM interfaces^ 

The theory of the anomalous T c (df rn ) dependence is 
based on the Usadel equations— describing the prox- 
imity effect of FM and SC layers, but neglecting a 
possible influence of the domains in the FM. In this 
manuscript we will show that an inhomogeneous stray 
field B stra y produced by the domain structure of a FM 
can actually also lead to a significant change in T c . To 
demonstrate this effect, we measure T c as a function of 
the perpendicularly applied magnetic field H of a Pb 
film on top of a Co/Pt multilayer with perpendicular 
magnetic anisotropy. In this sample the proximity effect 



is suppressed by an amorphous Ge layer between Pb and 
Co/Pt. The domain structure in the Co/Pt multilayer 
can consist of stable band or bubble domains. The FM 
layer can also be in a single domain state, depending on 
the preceding magnetization procedure, as was shown in 
a recent study of the vortex pinning in this system— i 
Due to field cancellation effects between H and B str ay, 
and due to the suppression of ip by B stra y, T C (H) can be 
controlled by changing the microscopic domain structure. 



II. MAGNETIC PROPERTIES OF THE Co/Pt 
MULTILAYER 

The properties of the Co/Pt multilayer have been 
described before— i Briefly, the multilayer has a 
[Co(0.4 nm)/Pt(1.0 nm)]i structure on a 2.8 nm Pt 
base layer on a Si/Si02 substrate. The magnetic proper- 
ties were characterized by the magneto-optical Kerr ef- 
fect (MOKE) and magnetic force microscopy (MFM), 
revealing that the sample has perpendicular magnetic 
anisotropy. Fig. ^ shows the magnetization Mf m of 
the Co/Pt multilayer measured by MOKE, normalized 
to the saturation magnetization M sat , as a function of 
the magnetic field H applied perpendicular to the sur- 
face. The loop has an almost rectangular shape with 
HoH n = 60 mT, hqH c — 93 mT, and fioH s = 145 mT, 
where H n , H c and H s are the nucleation, coercive and 
saturation field, respectively, and /io is the permeability 
of the vacuum. 

Using different magnetization procedures, one can pro- 
duce different stable domain patterns in the sample. For 
instance, after out-of-plane demagnetization, band do- 
mains are observed by MFM, see Fig.^b). Stable bubble 
domains with local magnetic moments m either pointing 
up (m z > 0) or down (m z < 0) perpendicular to the sam- 
ple surface can be created by applying a negative field of 
— IT, sweeping H to a positive value between H n and 
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FIG. 1: Magnetic properties of the Co/Pt multilayer: (a) Hys- 
teresis loop measured by magneto-optical Kerr effect with 
H perpendicular to the sample surface. MFM images 
(5x5 ^m 2 ) show that the domain structure of the sample 
consists of band domains after out-of-plane demagnetization 
(b), bubble domains in the s = 0.3 (c) and s = 0.93 (d) states. 

H s , and then removing H. The parameter s, which gives 
the fraction of magnetic moments that are pointing up 
(m z > 0) to the total amount of magnetic moments, is 
used to describe the different remanent magnetic states 
obtained after this magnetization procedure. The value 
of s can be found from the MFM images by dividing the 
dark area (m z > 0) by the total area, or from measure- 
ments of Mf m as will be described later. 
The lateral size of the domain structures can be esti- 
mated from the MFM images. The typical diameter of 
the bubble domains is about ~ 300 nm. The same value 
is obtained for the average width of the band domains. 
Although the magnetic moments of the Co/Pt multilayer 
are equally distributed between up- and down-directions 
in both the demagnetized and the s = 0.5 state, there 
are distinct differences between these two domain states: 



In the demagnetized state the lateral size of the domain 
is larger (because band domains are extended in one di- 
rection). Note also that in the demagnetized state, the 
boundary between domains with magnetization pointing 
up and down is well defined, but not straight: several 
sharp corners of different angles can be seen. No MFM 
images could be obtained for the s = 0.5 state, caused 
by the difficult magnetization procedure due to the steep 
slope of the Mf m (H) curve, see Fig.^a). However, from 
the image of the s = 0.3 state, see Fig. ^c), one can ob- 
serve that the domain walls are less sharp defined than 
in the demagnetized state. 

III. PHASE BOUNDARY OF THE 
SUPERCONDUCTING FILM 

After characterizing the properties of the FM, a 10 nm 
Ge film, a 50 nm Pb film and a 30 nm Ge capping layer 
are subsequently evaporated on the Co/Pt multilayer at 
a substrate temperature of 77 K. The amorphous Ge film 
between Pb and Co/Pt is insulating at low temperatures, 
so that the proximity effects between Pb and Co/Pt are 
suppressed. 

The upper critical field H C 2 of bulk type-II SCs is given 
bj^ 

with $o = 2.068 mT /an 2 the superconducting flux quan- 
tum, £(T) = £(0)/\/l — T/T c0 the temperature depen- 
dent coherence length in the dirty limit, and T c0 the crit- 
ical temperature at zero field. Hence, the linear slope of 
H c 2 as a function of temperature is only determined by 
the coherence length £. 

H C 2(T) can behave differently when the geometry of 
the SC is changed, e.g. for thin films with thickness 
w < £(T). While eq. Q]rs still valid for thin type-II su- 
perconducting films with H applied perpendicular to the 
sample surface, T c for parallel H is given by«2& 

Tm=T c »[l- * 2 ^\ lH\ (2) 

with T c q the zero-field critical temperature. In fact, this 
formula also gives the phase boundary of a mesoscopic 
line in perpendicular field, because the cross section, ex- 
posed to the applied field, is the same for a film of thick- 
ness w in parallel H and for a mesoscopic line in per- 
pendicular H. 27 For multiply connected mesoscopic lines, 
T C (H) can show an even more complex behavior due to 
fluxoid quantization effects^ 

The phase boundary of the SC/FM bilayer was mea- 
sured in a Quantum Design superconducting quantum 
interference device (SQUID) magnetometer with H ap- 
plied perpendicular to the surface. Fig.[3]shows the data 
obtained in two field cooled measurements of the total 
magnetization M = Mf m + M sc (M sc is the magnetiza- 
tion of the SC) at the applied field of /j,qH = 0.5 mT, 
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FIG. 2: Field cooled M(T) measurements of the hybrid 
SC/FM bilayer in fi H = 0.5 mT with the Co/Pt multilayer 
in the s = 0.5 (•) and s = (o) states. The arrows indicate 
T c . 

after the samples were brought in the s = 0.5 and s = 
states. Above T c , M has a constant value for both states, 
given by the contribution of the FM Mf m , from which 
s can be derived. When the sample is cooled through 
T c , a diamagnetic response of the SC appears. These 
kinds of measurements were used to determine T C (H) as 
the temperature where M starts to deviate from M/ m . 
Repeating these measurements at several applied fields 
\H\ < 25 mT did not change the offset M/ m above T c , 
implying an unchanged domain state. 

A. T C (H) with magnetized FM 

The phase boundary for the s = state (all m z < 0) 
obtained by several M(T) measurements in varying fields 
is shown in Fig. Ola). A linear behavior of the phase 
boundary is observed, which can be fitted by eq. ^with 
£(Q) = (41.2 ±0.2) nm and T c0 = (7.227 ± 0.002) K. This 
implies that in this state, the FM has no influence on 
the superconducting film, because both the linear behav- 
ior and the values of T c0 and £ (0) are in good agreement 
with those of pure Pb films S It is important to note that 
the temperature dependence of £(T) = £(0)/yl — T/T c0 
derived for this domain state is the same for all domain 
states, since we are always dealing with the same Pb film. 
Let us consider the magnetic stray field B stray of a homo- 
geneously magnetized film in the s — state, schemati- 
cally drawn in Fig. 01b) . B stra y has its largest amplitude 
at the sample boundary and is negligible above the center 
of the FM. Intuitively this can be understood by consid- 
ering the stray field of a single magnetic dipole in the 
center of the sample: The negative field above the dipole 
is compensated by the returning positive stray field of the 
surrounding magnetic dipoles. Therefore, the main cen- 
tral part of the superconductor is only weakly influenced 
by B s tray, and the measured T C (H) curve resembles the 




FIG. 3: (a) Magnetic field - temperature phase diagram of 
the superconducting Pb film covering the Co/Pt multilayer. 
Before measuring T C (H), the Co/Pt multilayer was brought 
into the s — state, (b) Schematic drawing of the stray field 
B s tray when the FM is in the s = state. 



one of a single Pb film. 



B. T C (H) with demagnetized FM 

The phase boundary for the demagnetized state, cor- 
responding to the MFM image shown in Fig. nib), is 
shown in Fig. Ufa). In this state, T c0 is suppressed 
to (7.048 ± 0.002) K. Moreover, the phase boundary 
still shows a linear behavior, but with a slope increased 
by a factor of 2.4. The difference between the phase 
boundaries in the demagnetized and the s = state 
can be attributed to the influence of the stray field 
B s tray, suppressing the order parameter in the super- 
conductor above T = 7.048 K. The coherence length at 
this temperature is £ = 260 nm. This means that su- 
perconductivity nucleates when the value of £ becomes 
smaller than approximately the width of the band do- 
mains. The nucleation first takes place in regions of 
the Pb film where the effective field in the z-direction 
fioH e ff z = /J>oH z + Bstray. z is minimum. The confine- 
ment of these superconducting nuclei leads to the differ- 
ent T C (H) dependence compared to the magnetized state. 
Aladyshkin et al. have very recently calculated the T C (H) 
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FIG. 4: (a) Magnetic field - temperature phase diagram of the 
superconducting Pb film covering the Co/Pt multilayer. The 
Co/Pt multilayer was demagnetized before measuring T C (H). 
As a reference we have added the phase boundary of the s = 
state. The dashed lines are guides to the eye. (b) Schematic 
drawing of the stray field B s tray when the FM has been de- 
magnetized. 
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FIG. 5: (a) Magnetic field - temperature phase diagrams of 
the superconducting Pb film covering the Co/Pt multilayer. 
Before measuring T C (H), the Co/Pt multilayer was brought 
into the s = 0.1 and s — 0.85 states. 



two dimensions. 



phase boundary of similar systems as the one that is ex- 
perimentally investigated here in the framework of the 
linearized Ginzburg-Landau equation. 30 They found, in 
agreement with our experimental result, that the upper 
critical field of a superconducting film can have very un- 
usual temperature dependencies when a single domain 
wall or periodic domain structures are present in a ferro- 
magnetic film that is in contact with the superconductor. 
Based on these considerations, we can conclude that the 
increased slope of the T C (H) curve may be related to the 
specific domain pattern in the demagnetized state. 
The MFM image in Fig. ^b) shows an equal contrast 
above all bright or dark domains, indicating that B stra y 
is rather homogeneous above the domains, and inhomo- 
geneous above the domain walls, which define some sharp 
corners. From calculations of the upper critical field 
of mesoscopic superconducting structures, e.g., triangles 
or squares, it is well known that the nucleation of su- 
perconductivity takes place first in the corners of these 
structures^i* 3 ^ Therefore, when trying to calculate the 
T C (H) phase boundary in order to explain the increased 
slope, one should take into account these corners formed 
by the domain walls. This could be done by expanding 
the one-dimensional model used by Aladyshkin et to 



C. T C (H) with bubbles in the FM 

The phase boundary shown in Fig. Ufa) is obtained 
when the Co/Pt multilayer contains bubble domains. 
Fig. Eland Fig. 0] are symmetric with respect to H, i.e., 
T c is the same for positive or negative H, but the pres- 
ence of the bubble domains causes an asymmetry of T c 
with respect to H. For bubbles having positive magnetic 
moments, i.e. s < 0.5, a higher T c is observed for positive 
H than for corresponding negative H, whereas for bub- 
bles containing negative magnetic moments (s > 0.5), 
T c is higher for negative H. Moreover, both T C (H) 
curves shown in Fig.[5fa) show a non- linear behavior with 
bumps in the field ranges around \(J-oH\ ~ 5 — 10 mT. 
To explain the asymmetric T C (H) curves, let us assume 
that the sample contains bubble domains with m z > 
in a matrix of magnetic moments with m z < 0, as shown 
in Fig. GJb): B stra y,z is positive above the bubbles and 
negative between them. A positive H in the z-direction 
compensates the negative B stra y,z between the bubbles 
and enhances B stra y,z above them, while a negative H 
has the opposite effect: it enhances £>, 
bubbles and compensates B, 



stray ,z 



between the 



stray, z 



above them. The im- 



portant point that causes the asymmetric phase bound- 
ary is that the absolute value of B stra y.z is larger above 
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FIG. 6: Dependence of the critical temperature at zero field 
T c (H=0) on the parameter s. The minimum value of T c is 
observed for s = 0.5. 



the m z > regions (bubbles) compared to the m z < re- 
gions (between the bubbles). When the sample is cooled 
in positive H, superconductivity can nucleate at higher 
temperatures in the area between the bubbles (where 
B stray, z < 0), compared to cooling the sample in the 
corresponding negative H, where the nucleation takes 
place in the areas above the bubbles. Note that qualita- 
tively similar non-linear T C (H) curves as those presented 
in Fig. 03a) have also been predicted by Aladyshkin et 
alM 

The critical temperature of the superconductor decreases 
when the bubble domains have larger density. To illus- 
trate this effect, Fig.[S]shows the dependence of T C (H = 
0) on the parameter s. A clear minimum of T c is ob- 
served around s — 0.5, which indicates that this domain 
state has the largest value of the stray field of all inves- 
tigated domain structures. Note that T c of the s = 0.5 
state is even lower than T c of the demagnetized state, 
emphasizing the inhomogeneous character of the s = 0.5 
state. The phase boundary of this domain state will be 
discussed in the next section. 



FIG. 7: Magnetic field - temperature phase diagrams of the 
superconducting Pb film covering the Co/Pt multilayer. Be- 
fore measuring T C (H), the Co/Pt multilayer was brought into 
the s = 0.5 state. As a reference the phase boundary of the 
s = state is added. The dashed line is a guide to the eye. 



using eq. EJand £(0) = 41.2 nm (from the phase bound- 
ary for s — 0), we obtain values of w — (213 ± 6) nm 
and T c0 = (6.994 ± 0.003) K. The value of w deter- 
mined from this fit can be compared with the typical 
bubble domain size of ~ 300 nm. For [IqH > 15 mT 
and T < 6.90 K, T C (H) shows a crossover from the one- 
dimensional network like to a two-dimensional linear be- 
havior, because the assumption w < £(T) for eq. [21 is no 
longer fulfilled. This is in agreement with the value of 
£(6.90 K) = 194 nm. 



IV. CONCLUSION 



D. T C (H) with the FM in the s = 0.5 state 

The phase boundary of the SC with the FM in the 
s = 0.5 state is shown in Fig. [7| In this domain state, 
E 'stray has a more inhomogeneous character than in the 
demagnetized state. For a discussion of the differences 
between these two domain states we refer to section [H] 
T C (H) follows a non-linear behavior, in contrast to the de- 
magnetized and the s = states. T C (H) for the s = 0.5 
state can not be described by eq. ^ but rather by eq. El 
in fields hqH < 15 mT, see the fit in Fig. EJb). This 
indicates that in the s — 0.5 state the regions where su- 
perconductivity nucleates can be considered as supercon- 
ducting strips with a width w < £(T), forming a sort of a 
superconducting network. When fitting the T C {H) curve 



In conclusion, the phase boundary between the normal 
and the superconducting state of FM/SC bilayers has 
been found to be strongly dependent on the domain 
structure in the FM. The stray field B s tray of these 
domains can lead to a significant decrease of T c in zero 
applied field, but, on the other hand, it can also enhance 
T c in applied fields. It has been demonstrated that the 
presence of bubble domains leads to the formation of the 
field-polarity dependent asymmetric phase boundaries 
T C (H) with respect to H, due to compensation effects 
between H and B stray For a specific inhomogeneous 
domain structure, the T C (H) phase boundary shows a 
crossover from a one-dimensional to a two-dimensional 
nucleation behavior. 
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